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The transition helicity amplitudes from the proton ground state to the Pn(1440) and Di3(1520) 
excited states (-y v pN* electrocouplings) were determined from the analysis of nine independent one- 
fold differential tt + tt~p electroproduction cross sections off a proton target, taken with CLAS at 
photon virtualities 0.25 GeV 2 < Q 2 < 0.60 GeV 2 . The phenomenological reaction model was 
employed for separation of the resonant and non-resonant contributions to the final state. The 
Pn(1440) and Di3(1520) electrocouplings were obtained from the resonant amplitudes parametrized 
within the framework of a unitarized Breit-Wigner ansatz. They are consistent with results obtained 
in the previous CLAS analyses of the 7r + n and tv°p channels. The successful description of a 
large body of data in dominant meson-electroproduction channels off protons with the same "/ v pN* 
electrocouplings offers clear evidence for the reliable extraction of these fundamental quantities 
from meson-electroproduction data. This analysis also led to the determination of the long-awaited 
hadronic branching ratios for the £>i 3 (1520) decay into Att (24%-32%) and Np (8%-17%). 

PACS numbers: 11.55.Fv, 13.40.Gp, 13.60.Le, 14.20. Gk 



I. INTRODUCTION 

An extensive research program on nucleon resonance 
(N*) excitation is in progress using the CLAS detector 
in Hall-B at Jefferson Lab The studies of tran- 

sition helicity amplitudes from the proton ground state 
to its excited states (or "f v pN* photo-/electrocouplings) 
represent a key direction in the N* program with CLAS. 
Meson-electroproduction data off nucleons in the N* re- 
gion obtained with CLAS open up an opportunity to de- 
termine the Q 2 -evolution of j v NN* electrocouplings in 
a combined analysis of various meson-electroproduction 
channels. The Q 2 -evolution of j v NN* electrocouplings 
will allow us to pin down active degrees of freedom in the 
N* structure at various distance scales and to access non- 
pcrturbativc strong-interaction mechanisms that govern 
the excited nucleon state formation as bound systems of 
quarks and gluons. 

Theoretical and experimental studies of the electroex- 
citation of nucleon resonances have a long history. Along 
with the hadron masses and their partial decay widths, 
the information on the ^ v pN* electrocouplings played an 
important role in the development of the quark mod- 
els in their contemporary advanced relativistic version in 
light-front dynamics @~14|. The picture of the nucleon 
and its excited states, which seemed quite simply mod- 
eled with three relativistic constituent quarks, turned out 
to be more complex. Recently obtained electrocouplings 
for the Pu(1440), P>i 3 (1520), and S u (1535) states [l| 
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showed that those quark models, which were successful 
in describing the electrocouplings of these states for Q 2 
> 2.0 GeV 2 , failed to reproduce the results at smaller 
photon virtualities of Q 2 < 1.0 GeV 2 . Quark models 
with flavor-conserving quark interactions are unable to 
describe the small mass of the Pn(1440) state [l6| and 
the ordering of Pn(1440) and 5n(1535) resonances [l[. 
Moreover, models that treat the Pn(1440) structure as 
just three constituent quarks are unable to describe the 
large total Pn(1440) decay width of 300 MeV. These 
difficulties prompted a search for additional contributions 
to the N* structure. 

A general unitarity requirement imposes mcson-baryon 
dressing contributions to both resonance electromagnetic 
excitation and hadronic decay amplitudes. Studies of 
meson-baryon dressing contributions to "f v pN* electro- 
couplings and resonance hadronic decay amplitudes, car- 
ried out at the Excited Baryon Analysis Center (EBAC) 
at Jefferson Lab (l7l-[24|. have extended our insight into 
the spectrum and structure of excited nucleon states con- 
siderably. The contributions from meson-baryon dressing 
to the Pn(1440) and £>i3(1520) electrocouplings were de- 
termined from a global analysis of the world data on ttN 
scattering and 7r + n, ir°p electroproduction off protons 
within the framework of the EBAC dynamical coupled- 
channel approach (EBAC-DCC) [H. 

This analysis showed that the contributions from 
mcson-baryon dressing to the j v pN* electrocouplings are 
maximal at small Q 2 and decrease with increasing pho- 
ton virtualities [13 • At Q 2 < 1.0 GeV 2 these contribu- 
tions may even be dominant. The meson-baryon cloud 
has a profound impact on the resonance spectrum. For 
example, in the Pn partial wave a single bare resonance 
pole located at W=1.76 GeV, being affected by meson- 
baryon dressing, splits into three poles located on differ- 
ent Riemann sheets [2l], [22[. Two of them with Re(W) 
« 1.36 GeV correspond to the physical Roper resonance. 
The double-pole structure of the Roper resonance was 
also observed in previous studies [25j, [26| . Meson-baryon 
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dressing should be taken into consideration in the inter- 
pretation of resonance electrocouplings as well as for the 
excited nucleon spectrum. 

Our studies of resonance electrocouplings at small pho- 
ton virtualitics presented in this paper offer valuable in- 
formation to further explore the role of meson-baryon 
and quark components in the N* structure. The sepa- 
ration between the meson-baryon cloud and quark core 
contributions within a well-defined theoretical framework 
[27j can help pin down the domain of photon virtualitics 
where quark components are the main contributor to the 
N* structure. This kinematic domain is of particular 
interest for the studies of hadrons from the first princi- 
ples of QCD, including Lattice QCD (LQCD) and 
Dyson- Schwinger equation studies of QCD (DSEQCD) 

The CLAS detector at Jefferson Lab is a unique large- 
acceptance instrument designed for the comprehensive 
exploration of exclusive meson electroproduction. It of- 
fers excellent opportunities of studying the electroexcita- 
tion of nucleon resonances in detail and with precision. 
The CLAS detector has provided the dominant portion 
of all data on meson electroproduction in the resonance 
excitation region. 

A variety of measurements of single pion electro- 
production off protons, including polarization measure- 
ments, have been performed at CLAS in a range of Q 2 
from 0.16 to 6 GeV 2 (42|. The electroexcitation ampli- 
tudes for the low-lying resonances p33(1232), Pn(1440), 
Pi3(1520), and 5*11(1535) were determined over a wide 
range of Q 2 in a comprehensive analysis of JLab-CLAS 
data on differential cross sections, longitudinally po- 
larized beam asymmetries, and longitudinal target and 
beam-target asymmetries |l5|. 

The combination of the large-acceptance CLAS detec- 
tor and the continuous electron beam from CEBAF made 
it possible to measure tt + tt~p electroproduction cross sec- 
tions with nearly full kinematic coverage for this three- 
body final hadron state @, |43| . These are the most ex- 
tensive data sets on unpolarized n + ir~p electroproduc- 
tion cross sections obtained so far. These data allowed 
for the first time the projection of nine one-dimensional 
differential cross sections, each sensitive to a different 
combination of resonance and background strength. The 
data of 0] were collected in the mass range 1.31 GeV < 
W < 1.56 GeV and with photon virtualities 0.25 GeV 2 
< Q 2 < 0.6 GeV 2 . A good description of these data 
was achieved within the framework of a phenomeno- 
logical Jefferson Laboratory - Moscow State University 
( JM) reaction model Q , that allowed us to establish the 
mechanisms contributing to this exclusive reaction. The 
presence and strength of the contributing tt + tt~p elec- 
troproduction mechanisms were established by studying 
their kinematical dependencies and correlations in differ- 
ent one-fold differential cross sections. 

In this work we present results on the electroexcitation 
of the Pn(1440) and £>i 3 (1520) states, obtained from 
the analysis of data on tt + tt~p electroproduction off pro- 



tons @ . The analysis was carried out employing the JM 
reaction model 0, which was further developed to pro- 
vide a framework for the determination of ^ v pN* electro- 
couplings from a combined fit of unpolarized differential 
cross sections. In the previous studies Q we did not at- 
tempt to isolate the contributions from resonances. In 
the analysis reported in this paper, we employ the JM 
model with the goal of isolating the resonant contribu- 
tions for the individual differential cross sections. For 
the description of resonant amplitudes, we updated the 
Breit-Wigner (BW) parametrization, making it consis- 
tent with the restrictions required by the general uni- 
tarity condition. The reliable evaluation of the resonant 
contributions enabled for the first time the determina- 
tion of the j v pN* electrocouplings for the Pn(1440) and 
.Di3(1520) states from charged double pion electropro- 
duction off protons. This complements the results from 
Ntt electroproduction in an independent channel. 

Analyses of different exclusive channels are essen- 
tial for a reliable extraction of resonance parameters. 
Currently the separation of resonant and non-resonant 
parts of the electro production amplitudes can be done 
only within phcnomcnological reaction models. There- 
fore, the resonance parameters extracted from the me- 
son electro production data fit may be affected by the 
model assumptions, and their credibility should be fur- 
ther examined. Non-resonant mechanisms in various 
meson-electroproduction channels are completely differ- 
ent, while the *y v NN* electrocouplings are the same. In- 
dependent analyses of different exclusive channels make 
it possible to test whether they give consistent results 
for the resonance electrocouplings. Most nucleon reso- 
nances decay into both Ntt and Ntttt final states. Stud- 
ies of resonance electroexcitations in these channels with 
completely different non-resonant contributions offer in- 
dependent information on N* electrocouplings. There- 
fore, a successful description of the data on tt + ji, ir°p, 
and tt + tt~p electroproduction off protons with consis- 
tent N* electrocoupling values provides clear evidence 
for the reliable extraction of these quantities from meson- 
electroproduction data. 

Studies of the Pi 1 (1440) resonance in tt + tt~p electro- 
production off protons offer the additional opportunities 
to improve the knowledge on electrocouplings of this ~ 
300 MeV broad state. Contrary to Nir electroproduction 
channels, the P33(1232) resonance does not directly con- 
tribute to the resonant parts of the tt + tt~p electroproduc- 
tion amplitude. Hence, the influence of the pj3(1232) res- 
onance on the extracted electrocouplings of the Pn (1440) 
is much weaker. 

The requirement of (^-independent N* hadronic de- 
cay amplitudes in tt + tt~p electroproduction provides con- 
straints on the Ntttt partial decay widths. This makes 
possible access to the ttA and pp partial hadronic decay 
widths of the Pn(1440) and L»i 3 (1520) in the measure- 
ments of tt + tt~p electroproduction off protons. 
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II. PHENOMENOLOGICAL MODEL JM FOR 
EVALUATION OF y v pN* ELECTROCOUPLINGS 

The phenomenological meson-baryon model JM was 
developed to describe tt + it~p electroproduction off pro- 
tons 0, mil with the primary objective of determin- 
ing the resonance "/ v pN* electrocouplings and the 7rA 
and pp partial hadronic decay widths from a combined 
fit to all measured observablcs. In our current analy- 
sis of the CLAS tt + tt~p electroproduction data @ the 
JM model was used to separate the resonant and non- 
resonant contributions to differential cross sections and 
to access the electrocouplings and 7rA, pp decay widths of 
the Pn(1440) and Di3(1520) resonances. Here we briefly 
discuss the basic ingredients of the JM model that are 
relevant to the objectives of this paper. 



A. Kinematics and cross sections 

At a given invariant mass W and photon virtuality 
Q 2 , the 7„p — > ir + Tr~p reaction can be fully described as 



a 
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FIG. 1: Kinematic variables for the description of ep — > 
e'p'-K + ir~ in the CM frame of the final-state hadrons corre- 
sponding to the first assignment presented in Section III Al 
Panel (a) shows the n~ spherical angles 9 W ~ and ip w - . Panel 
(b) shows the angle a^-j^/^+j between the two planes: one 
of them (plane A) is defined by the 3-momenta of the ini- 
tial proton and the final tt~ , the other (plane B) is defined 
by the 3-momenta of the two others final hadrons tt + and 
proton. The unit vectors 7 and f3 are normal to the ir~ three- 
momentum in the planes A and B, respectively. 



a five-fold differential cross section d 5 a/d 5 r, where d 5 r is 
the phase-space volume of the five independent variables 
in the center-of-mass (CM) system of the final Tr + ir~p 
state. There are many possible choices [4£| of the five 
independent variables. After defining M^+ p , M w - p , and 
A/jt+tt- as invariant mass variables of the three possible 
two-particle pairs in the tt + it~p system, we adopt the 
following three assignments: 



1 



d 5 Ti = dM pn +dM v + 7r -dil v -dai p i 7r +^ 7r -^ where 
(#„.-, n — ) are the final tt~ spherical angles 
with respect to the direction of the virtual photon, 
and Q![p'7r+][p7r-] is the angle between Plane B de- 
fined by the momenta of the final p'7r + pair and 
Plane A defined by the momenta of the initial pro- 
ton and the final tt~ : 



2. d 5 r 2 = dM^+dMn 



dVL p idu\. 



][v'p) 



where fL 



(6 p > , tfpi ) are the final proton spherical angles with 
respect to the direction of the virtual photon, and 
a[ w + w -i[p/pi is the angle between Plane B defined 
by the momenta of the tt + it~ pair and Plane A de- 
fined by the momenta of the initial and final pro- 
tons : 



3. d 5 r 3 = dM p 



.+ dM p7T -dVL m 



-do,\. 



where f2„ 



(9^+ , tp v + ) are the final 7r + spherical angles with 
respect to the direction of virtual photon, and 
a [p'7i-][pTr+] is the angle between Plane B defined 
by the momenta of the final p'ir~ pair and Plane 
A defined by the momenta of the initial proton 
and the tt + . 

All frame-dependent variables are defined in the final 
hadron centcr-of-mass (CM) frame. 

The emission angles for the final state particles in the 
case of the first assignment are shown in Fig. [1] This 
choice is most suitable for describing Tt + n~p electropro- 
duction through the 7r~A ++ intermediate state, which 
is the dominant contributor of all isobar channels in the 
kinematic region covered by the data Q. For the other 
assignments the emission angles of the final hadrons are 
analogous to the ones given in Fig. [1] The relations be- 
tween the momenta of the final-state hadrons and the 
five variables of the first assignment can be found in the 
Ref. @. 

The 7r + 7r~p electroproduction data have been collected 
in the bins of a seven dimensional space. As mentioned 
above, five variables are needed to fully describe the fi- 
nal hadron kinematics, while to describe the initial state 
kinematics two others variables W and Q 2 are required. 
The huge number of seven dimensional bins over the re- 
action phase space (ps 100,000 bins) does not allow us 
to use the correlated multi-fold differential cross sections 
in the analysis of electroproduction processes, where the 
statistics decrease drastically with the photon virtualitics 
Q 2 . More than half of the five-dimensional phase-space 
bins of the final hadrons is not populated due to statis- 
tical limitations. This is a serious obstacle for any anal- 
ysis method that employs information on the behavior 
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FIG. 2: Reaction mechanisms of the JM model Q that contribute to 7r + 7r p electroproduction in the kinematic region covered 
by recent CLAS measurements [(| with W < 1.6 GeV and 0.25 GeV 2 < Q 2 < 0.6 GeV 2 : 3-body mechanisms (a), ttA isobar 
channels (b), and direct 2-7T production mechanisms that correspond to different assignments for the final state hadrons (c,d). 
The ti (i=l,..3) stand for the squared transferred momenta in the exchange processes by unspecified particle(s), as described 
in Q and shown by blobs in the (c) and (d) panels. 



of multi-fold differential cross sections. We therefore use 
the following one-fold differential cross sections in each 
bin of W and Q 2 covered by the measurements: 

• invariant mass distributions for three pairs of 
the final particles da / dM n +^- , da/dM^+„, and 



da/dM^-p] 

• angular distributions for spherical angles 
of the three final particles da/d(—cos0^-), 
do / d(—cos9 7r +)- l and da / d{—cos6 p i) in the CM 
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frame; 

• angular distributions for the three a-anglcs 
described above and determined in the CM 
frame da/da^^+^p^-], da/da[ p i 7T -]ipn+] > an d 
da/da^+n-Wpp'}- 

The one-fold differential cross sections were obtained 
from integrating the five-fold differential cross sections 
over the relevant four others kinematic variables of d 5 Ti. 
All details related to the evaluation of the tt + tt~p one- 
fold differential cross sections we are using for the extrac- 
tion of resonance parameters can be found in Ref. . 

B. Relevant electroproduction mechanisms 

The major part of the ir + n~p electroproduction off 
protons at W < 1.6 GeV is due to contributions from 
the two irA isobar channels ir~A ++ and 7r + A°. The 
A ++ (1232) resonance is clearly seen in all ir + p mass dis- 
tributions for W > 1.4 GeV, while contributions from 
the 7T + A° isobar channel arc needed to better describe 
the data in the low mass regions of the w~p mass dis- 
tributions. The observed in the data [f| strength of 
the 7r~A ++ isobar channel is approximately nine times 
larger than that of 7r + A° due to isospin invariance. 
The contributions from all other isobar channels p p, 
7r+L>? 3 (1520), 7r+^ 1 ° 5 (1685), and tt" K^+(1640), which 
arc incorporated into the JM-model [45l448| in order to 
describe the data at W > 1.6 GeV, are negligible in the 
kinematic region covered in this analysis, and are not 
included in this work. 

The production amplitudes relevant for our analysis 
of ^ v p — > ir + ir~p are illustrated in Fig. [2] They con- 
sist of the 7T~A ++ and 7r + A° isobar channels and di- 
rect double pion production mechanisms. The produc- 
tion amplitudes for irA intermediate states (Fig. [2b) con- 
sist of the resonant contributions -f v N — > N*, A* — > 7rA 
and non-resonant terms. All resonances listed in Ta- 
ble |T] are included in the JM model. However, in the 
kinematic area covered in our measurements, only the 
Pn(1440), Di3(1520), and Sii(1535) nucleon resonances 
have strength that is sufficient to manifest themselves 
in the one-fold differential cross-sections. Non-resonant 
amplitudes, depicted in the diagrams in Fig. [2b are com- 
puted from the well-established Born terms presented in 
Appendix A of Ref. 0. The additional contact terms 
are implemented in the full ttA production amplitudes. 
They describe effectively the contributions from mecha- 
nisms other than the Born terms to it A production and 
the part of the irA final-state interactions (FSIs) that 
are not included in the JM model's absorptive approxi- 
mation for FSIs 0, H3] ■ The Lorcntz structure of these 
additional contact terms is determined by superposition 
of the two second-rank Lorentz tensors 

P S ds9^, (1) 
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-Di 3 (1700) 


1.74 
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TABLE I: List of resonances invoked in the tt + tv~p fit and 
their parameters: total decay widths T to t and branching frac- 
tions (BF) to 7rA and pp final states. The quoted values 
for the hadronic parameters are taken from fits to the ear- 
lier CLAS -k + ty~p data [|3] using the 2005 version of the JM 
model [46l . [4?J . The quantities labeled as var correspond to 
the variable parameters fit to the CLAS 7r + 7r~p data 0] within 
the framework of the current JM model version 0] employ- 
ing the unitarized BW ansatz of Section Til Dl for the resonant 
contributions. Start values for the resonance electrocouplings 
are taken from the references listed in the last column and 
extrapolated to the Q 2 area covered by the CLAS experiment 



where p c = (2p™ — q-y) is the difference of the final pion 
four-momentum pi and the momentum transferred, q 7 — 
p^,. A parametrization of the additional contact-term 
amplitudes in the tt A channels may be found in Appendix 
B of Ref. 0. 

All isobar channels combined account for 70% to 90% 
of the charged doublc-pion fully integrated cross sections 
in the kinematic region covered by data @ . The remain- 
ing part of the cross sections comes from direct charged 
pion (2ir) production mechanisms, in which the tt + tt~p 
final state is created without the formation of unstable 
hadrons in the intermediate states. The presence of these 
mechanisms is required by the unitarity of the three- 
body ir + iT~p production amplitudes [Edl ]. Their manifes- 
tation in n + iT~p electroproduction was observed for the 
first time in our previous analyses of CLAS data 0, [47j . 
The dynamics of these processes was unknown and has 
been established from the CLAS data analysis within the 
framework of the JM model. The direct 27r-production 
mechanisms incorporated in the JM model are depicted 
in Fig. [2Jc,d). They represent two subsequent exchanges 
of unspecified particles, parametrized by propagators 
that depend exponentially on the running four-momenta 
squared. Each set of diagrams in Fig. [2] corresponds to 
various assignments of the final-state hadrons, resulting 
in different four-momenta squared running over propa- 
gators in the exchange amplitudes. The JM model ex- 
tends for the first time the description of tt + tt~p electro- 
production beyond the approximation of superimposed 
isobar channels by incorporating direct 27r-production. 
Explicit expressions for the above-mentioned direct 27T- 



7 



production amplitudes can be found in Appendices A-C 
of Ref. 0- 

The relationships between Tr + ir~p clectroproduction 
cross sections and the three-body production amplitudes 
employed in the JM model are given in Appendix D of 
Ref. [tJ • This information is required in order to compare 
the amplitudes of the JM model with the results from any 
other study of ir + ir~p electroproduction amplitudes. 



C. Breit-Wigner parametrization of resonant 
amplitudes 

We start from a non-unitarized relativistic Breit- 
Wigner (BW) ansatz to describe the resonant contribu- 
tion (Xf \T r \ A 7 A P ) in the helicity representation: 



(A / |T r |A 7 A p )=^ 



(A/ \Tdec\ A_r)(A_r |T em | A 7 A P ) 
M 2 - W 2 - iT r (W)M r ' 



(2) 



where M r and T r are the resonance mass and energy- 
dependent total width, respectively. The matrix elements 
(A_r \T em \ A 7 A P ) and (A/ \T dec \ X R ) arc the electromag- 
netic production and hadronic decay amplitudes of the 
N* with helicity X R = A 7 — A p , in which A 7 and A p stand 
for the helicities of the photon and proton in the ini- 
tial state, and A/ represents the helicities of final-state 
hadrons in the N* decays. 

The hadronic decay amplitudes (Xf\Tdec\Xfi) are re- 
lated to the T\ f (W) partial hadronic decay widths of 
the N* to 7r A or pp final states / of helicity A / by: 

(X f \T dec \X R ) - (\ f \Tt\\ R Kl(cose*)e^\ 

with n = Xr and v = — A a for the 7rA intermediate state 
and v = X p > — X p for the pp intermediate state, and 



The means and (pi) are the magnitudes of the three- 
momenta of the final ir for the N* — > 7rA decay or 
of the final proton for the TV* — > pp decay (i=2), evalu- 
ated at W = M r and at the running W, respectively, and 
averaged over the running mass of the unstable hadron 
in the intermediate state. The variables 8*, </>* are the 
CM polar and azimuthal angles for the final ir (N* — > 7rA 
decay) or the final proton (N* — > pp decay). The sym- 
bol J r stands for the N* spin. This relationship between 
the N* hadronic decay amplitudes and the partial decay 
widths was derived for the N* — > ttA decay in our previ- 
ous article • All details on the parametrization of the 
resonance hadronic decay amplitudes in the JM model 
can be found in Appendix A of this paper. 

Energy dependencies of the partial hadronic decay 
widths are described under the assumption that centrifu- 
gal barrier-penetration factors are the major contributors 




to the off-shell behavior of resonance hadronic decay am- 
plitudes dillil: 



LS 



M r (J 2 {p r R) + N 2 {p r R)) 
W (J 2 (pR)+N 2 (pR)) 



1/2 



(4) 



where <// and Ni are the Bessel and Neumann func- 
tions. The factor in square brackets represent the ratio 
of barrier penetration factors for the final meson with 
orbital angular momentum I evaluated at W = M r and 
at the running W. The variable R represents the in- 
teraction radius whose value was set to 1 fm, and y / T£s 
stands for the decay amplitude estimated at the resonant 
point. The partial decay amplitudes VTls of Eq. ((4]) arc 
transformed from the LS to the helicity representation 
(see the Appendix A) and used to compute the ./V* 
hadronic decay amplitudes via Eq. 

In the JM model, the N* total decay width T r (W) 
in Eq. © is evaluated as a sum over all partial decay 
widths. In this way we ensure unitarity of the resonant 
amplitude for a single resonance contribution. The unita- 
rization procedure in the actual case of many contribut- 
ing resonances will be discussed in Section III PL 

The resonance electroexcitation amplitudes 
(Ah \T em \ A 7 A P ) in Eq. © are related to the -y v NN* 
electrocouplings A1/2, ^3/27 and S'1/2 for nucleons. The 
definition of these electrocouplings in the JM model is 
consistent with the Review of Particle Physics (RPP) 
[53l ] relation between the Ai/ 2 , A 3 / 2 electrocouplings 



and the N* electromagnetic decay width T n 



r 7 = 



17, r 



2M 



N 



(2J r + l)M r 



\A 1/2 



A 



3/2 1 



(5) 



where q 1 r is the three-momentum modulus of the photon 
at W — M r in the CM frame. The transition amplitudes 
(Xr \T em \ A 7 Ap) are related to the j v NN* Ai/ 2 , A 3 / 2 , 
and S'1/2 electrocouplings by imposing the requirement 
that the BW parametrization 54| of the resonant cross 
section for a single contributing resonance should be re- 
produced: 



a{W) 



^( 9 j 1 n M^mr, g 7 
q* { r 1 (M r 2 - W 2 ) 2 - M 2 T 2 (W) R- 



(0) 



Here the photon equivalent energy K = — 2 Jf N , q-y is 
the absolute value of the initial photon three-momentum 

of virtuality Q 2 > g 7 = yjQ 2 + E 2 and the energy _E 7 

in the CM frame at the running W 



E-y 



W 2 -Q 2 



M 2 



2W 



(7) 



The qy^ r value in Eq. §5§ can be computed from Eq. 
([7]) at W=M r . Ti(W) is the energy- dependent hadronic 
decay width to the final state 7rA (i=l) or pp (i=2). 
The factor ^ in Eq. © is equal to unity at the pho- 
ton point. It accounts for the choice @ of the virtual 
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photon flux in the evaluation of the virtual photon cross 
sections. In this way we obtain the following relationship 
between the transition amplitudes (Xr \T em \ A 7 A P ) and 
the ry v NN* electrocouplings: 

(\ R \T em \x,x p )-— y 4 ^ ^—a 1/2 . 3/2 (q ), 

1 3 

with |A 7 — Xp \ = — , — for transverse photons, and 

<A*|T ero |A 7 AA__y — sj — Si/ziQ ), 

for longitudinal photons. 

(8) 

The factor Ana in Eqs. © reflects the particular rela- 
tionship between the JM model amplitudes and cross sec- 
tions p} , when the absolute value of the electron charge 
is factorized out of the production amplitudes. 



D. Unitarization of the full resonant amplitudes 

The BW ansatz, described in Section III C[ provides 
unitary contributions from any individual TV* state to 
the resonant part of the electroproduction amplitudes. 
However, the full resonant amplitude, which represents 
a superposition of contributions from all TV* states, is 
not unitary. Figure [3] schematizes the processes that con- 
tribute to the resonant propagator dressing. They consist 
of transitions between the same and different TV* states 
mediated by the strong interaction. The regular BW 
ansatz incorporates only transitions between the same 
TV* states shown by the blob in the top panel of Fig. [3J 
The unitarization is achieved via inclusion of both dia- 
grams in Fig. 3, allowing us to take into account tran- 
sitions between both the same and different TV* states 
(bottom panel in Fig. [3]) in the dressed resonant prop- 
agators. We replaced the usual BW ansatz by its uni- 
tarized extension originally proposed in Ref. [551 ] ■ This 
ansatz was adjusted for the description of the BW reso- 
nance electroproduction amplitudes of Eq. (f2J) employed 
in the JM model. The full resonant amplitude of the uni- 
tarized extension of the BW ansatz satisfies the unitarity 
condition [55J. 

Unitarized resonance amplitudes in the JM model in- 
corporate photo- and electro-excitation amplitudes of all 
relevant resonances a, all possible transitions between 
the initial TV* state a and the final TV* state /3, and 
hadronic decays of j3. Fully unitarized resonance am- 
plitudes are determined by the sum of the products of 
the electromagnetic excitation amplitudes of the a-th TV* 
state {X a \T em \ A 7 A P ), the hadronic decay amplitudes of 
the (3-th TV* state (A/ \Td ec \ Xp), the dressed propagator 



M 




FIG. 3: Diagrams showing the unitarized BW ansatz of the 
JM model, which incorporates transitions via the dressed res- 
onance propagator between states listed in Table U for the 
same N a =Np (top) and different N a ^ Np resonances as al- 
lowed by the conservation laws of the strong interaction (bot- 
tom). 

S a p, and 

(A/ |T res | A 7 A P ) = ^](A/ \Tdec\ Xp)S a /3(X a \T em \ A 7 A P ). 

a,/3 

(9) 

The sum incorporates all transitions between the a-th 
initial and the /3-th final TV* states that are allowed by 
quantum number conservation laws. It runs over all re- 
peated indices a and /3 that label the TV* states. As a con- 
sequence of angular momentum conservation A a =A^, but 
these two helicities stand either for the same TV* states 
if a=/3 (diagonal transitions) or different TV* states if 
a 7^ (3 (off-diagonal transitions). The index / represents 
the final-state helicities, either for irA or for pp. 

The expression for the inverse dressed propagator S~p 
is obtained in Ref. [H[ (see Eq. (6) in Ref. HH). However, 
the parametrization of the TV* propagator in Ref. [55| 
and in JM are different. Therefore, we change the in- 
verse resonant propagator S~l of Ref. [Hj], so that it 
coincides with the single TV* contribution employed in 
the JM model. In this way we obtain 

S~l = M 2 N ,5ap - i(£ y/rU^)yjM N .M N ^ (10) 

k 

-W 2 5 a0 , 

where the index k represents the partial TV* hadronic 
decay widths to all possible final states, decomposed over 
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LS partial waves. For a single TV* contribution, Eq. (fTU)) 
coincides with the inverse BW propagator in Eq. ([2]). 

The resonant contribution at W < 1.6 GeV incorpo- 
rates the N* states listed in Table Hand the 3/2+ (1720) 
candidate state observed in the previous analysis of the 
CLAS tt + it~p electroproduction data [43j]. The conser- 
vation laws in strong interactions allow transitions in 
dressed resonant propagators only between the following 
pairs of N* states: D 13 (1520) and D 13 (1700), 511(1535) 
and 5u(1650), and 3/2+(1720) candidate and P 13 (1720) 
(bottom panel in Fig. [3]) . Therefore, the unitarized reso- 
nant amplitude represents the sum of regular BW ampli- 
tudes over all other TV* states with only diagonal tran- 
sitions in the dressed resonance propagators plus contri- 
butions from the aforementioned pairs of N* states that 
have both diagonal and non-diagonal transitions in the 
resonant propagators. The contributions from each pair 
of TV* states (A/ lT2r.es | A 7 A P ) are determined by Eq. © 
with indices a and /3 running from 1 to 2, resulting in 
2x2 S~l matrices. After invertion of S~p in Eq. (JTUJ) and 
insertion into Eq. ©, we get the following expressions 
for the contribution from each of the three pairs of N* 
states listed in Fig. [3J 

(A/ \T 2r es\ A 7 A p ) = 

det[S ap ] 

{(A/ \T dec \ Ai)(M| - iT 2 (W)M 2 - VF 2 )(Ai |T em | A 7 A P ) 
+ (A/ \T dec \ A 2 )(M 2 - iT 1 (W)M 1 - W 2 )(X 2 \T em \ A 7 A P ) 

+i(X f \T dec \ A 2 ) J2 Vri t r 2fc y/M N ;M N ;( Ai |T em | A 7 A P ) 

k 

+i(\ f \T dec \ Ai) £ y/T^yjM N ;M NS (\ 2 \T em \ A 7 A P )}. 
k 

(11) 

The first two terms in Eq. (fTTj) correspond to the diagonal 
S a p elements. They describe the processes depicted in 
the top panel of Fig. [3J The two other terms correspond 
to transition between different N* states, shown in the 
bottom panel of Fig. [3J The determinant of S~l can be 
computed from Eq. (fTUj) as 

det[S~l] = (AI 2 - iri(W)Afi - w 2 y 
(AI 2 - iY 2 (W)M 2 - W 2 ) + (J2 V^TT^AhAh. 

(12) 

It follows from Eqs. (jlip and (fT2"]) . that if transitions 
between two different N* states in the dressed resonant 
propagators become impossible, then 

1^ = 0- (13) 

k 

namely, all off-diagonal terms become zero, while the di- 
agonal terms give rise to the amplitudes of the regular 
BW ansatz. When the transitions between the pairs of 



N* states are turned on, they not only add off-diagonal 
terms to Eq. (fTTj) , but they also change the effective prop- 
agator in the (A/ |T2 res | A 7 A P ) amplitude with respect to 
the BW ansatz Eq. ©, determined by l/det[S~l]. 



III. THE CLAS DATA FIT 

The resonance parameters obtained in this paper are fit 
to the CLAS tt + tt~p electroproduction differential cross 
sections @. A realistic evaluation of resonance electro- 
coupling uncertainties extracted in the fit is an important 
objective of our analysis. We require that the ranges 
of resonance electrocouplings extracted from the 7r + 7r~p 
electroproduction channel take into account both uncer- 
tainties in the measured differential cross sections and in 
the JM model parameters. In order to provide a real- 
istic evaluation of resonance parameters, we abandoned 
the traditional least-squares fit, since the parameters ex- 
tracted in such a fit correspond to a single presumed 
global minimum, while the experimental data descrip- 
tion achieved with others local minima may be equally 
good within the data error bars. Furthermore, the tra- 
ditional evaluation of fit-parameter uncertainties, based 
on the error propagation matrix, cannot be used for the 
same reason. 

A special procedure was developed to obtain not only 
the best fit, but also to establish bands of computed cross 
sections that are compatible with the data within their 
uncertainties. In the fit we vary simultaneously non- 
resonant and resonant parameters of the JM model given 
in Tables U and [TT1 respectively, around their start val- 
ues, employing unrestricted normal distributions. The 
choice of the start parameters and the normal distribu- 
tion cr-values employed for parameter variation will be 
further described below. For each trial set of the JM- 
model resonant and non-resonant parameters we com- 
pute nine measured one- fold differential T: + Tr~p cross sec- 
tions, and the x 2 per data point values (x 2 /d.p.). The 
X 2 /d.p. values were estimated in point- by-point compar- 
isons between the measured and computed one-fold dif- 
ferential cross sections in all bins of W and Q 2 covered by 
the CLAS tt + tt~p data of Ref. @. In the fit we selected 
computed one-fold differential cross sections closest to 
the data with x 2 /d.p. less than a predetermined maxi- 
mum value Xmax./ d -P-- The values of x 2 n ax./ d -P- were 
obtained by requiring that the computed cross sections 
with smaller x 2 /d.p. be within the data uncertainties for 
the majority of the data points, based on point-by-point 
comparisons between the measured and the computed 
cross sections (see examples in Figs.lHE]). I n this fit pro- 
cedure we obtain the x 2 1 'd.p. intervals within which the 
computed cross sections describe the data equally well 
within the data uncertainties. 

The resonance parameters obtained from all these 
equally good fits are averaged and their mean values are 
taken as the resonance parameters extracted from the 
data. Dispersions in these parameters are taken as the 



10 





Ranges covered 




in variations of the 


Variable parameters 


start parameters, 
% from their values 


Magnitude of the additional 




contact term amplitude in the 


45.0 


7T~ A ++ sub-channel 




Magnitude of the additional 




contact term amplitude in the 


60.0 


7T + A° sub-channel 




Magnitudes of six 




2tv direct production 


30.0 


amplitudes 





TABLE II: Variable parameters of non-resonant mechanisms 
incorporated into the JM model 0. The ranges in the table 
correspond to the 3cr areas around the start values of param- 
eters. 



N* states 


Mass, 
(MeV) 


Total decay width, r to t, 
(MeV) 


Pn(1440) 
D 13 (1520) 
5n(1535) 


1430-1480 
1515-1530 
1510-1560 


200-450 
100-150 
100-200 



TABLE III: Allowed N* hadronic parameter variations in the 
fit of the CLAS -k + tv~p electroproduction data [(J. 



uncertainties. Our fitting procedure allows us to obtain 
more realistic uncertainties of the fit parameters than 
from the usual least-squares method. In this way we take 
into account both statistical uncertainties in the data and 
systematic uncertainties imposed by the use of the JM 
reaction model. 

We vary the parameters of the JM model, listed in 
Table ILT1 for the non-resonant mechanisms, which are rel- 
evant in describing the CLAS it + tt~p data @. They are 
limited to non-resonant irA sub-channels and direct 2w- 
electroproduction. 

The magnitudes of the additional contact terms in the 
7r~A ++ and ir + A° isobar channels (row (b) in Fig. [2]), 
as well as the magnitudes of all direct 27r-production am- 
plitudes (rows (c),(d) in Fig. [3]), are chosen as variable 
parameters. All these parameters were determined in the 
fit to the CLAS data on ir + ir~p electroproduction in our 
previous analysis without variation of the resonance 
parameters. In the CLAS data fit presented in this paper, 
the values of the aforementioned non-resonant parame- 
ters of the JM model are re-evaluated under simultaneous 
variation of: 

• the magnitudes of additional contact-term ampli- 
tudes in the 7r~A ++ and tt + A° isobar channels (2 
parameters per Q 2 bin); 

• the magnitudes of all direct 27r-production ampli- 
tudes (6 parameters per Q 2 bin); 

• the resonant parameters listed in Table HI The 
CLAS ir + ir~p data [1] are mostly sensitive to the 



electrocouplings of the Pu(1440) and £>i 3 (1520) 
states (5 resonance electrocouplings per Q 2 bin), 
to the 7tA and pp hadronic decay widths of these 
two resonances, and to the S , n(1535) state (6 pa- 
rameters that remain the same in the entire Q 2 area 
covered by the measurements). 

Therefore, we consistently account for the correlations 
between variations of the non-resonant and the resonant 
contributions while extracting the resonant parameters. 

The VF-dependences of the magnitudes of the addi- 
tional contact-term amplitudes in the 7rA sub-channels 
and of the magnitudes of the direct 2-7r-production am- 
plitudes are adopted from our previous analysis 0. 
We apply multiplicative factors to the magnitudes of 
the extra contact-term amplitudes and the 27r direct- 
production amplitudes. The multiplicative factors arc 
W-independent within any Q 2 -bin, but they are fit to 
the data in each Q 2 -bin independently. In this way we 
retain the smooth VF-dependences of the non-resonant 
contributions established in our previous analysis 0. 

We use two parameters for the variation of the mag- 
nitudes of additional contact-term amplitudes in the 
7r~A ++ and 7r + A° isobar channels. However, the non- 
resonant parameters in the 7r + A° isobar channel have a 
rather small impact on the fit, since the non-resonant 
contributions of the 7r + A° isobar channel are approxi- 
mately a factor 9 smaller than that of the 7r~A ++ . 

The parameters for the non-resonant Born terms in the 
7rA sub-channels ((b)-row in Fig. [5]) include the pirA cou- 
pling and cut-off for this hadron transition form factor, 
as well as the electromagnetic pion and nucleon form fac- 
tors (Hoj . All these parameters arc taken from previous 
studies of meson photo-, electro-, and hadroproduction 
referred to in Ref. and are kept fixed for the current 
fit. 

Different assignments of final hadrons in 27r-direct pro- 
duction mechanisms, as shown in Fig. [5] (rows (c),(d)), 
result in 12 different sub-processes. Our previous anal- 
ysis demonstrated that the magnitudes of all these 
sub-processes are determined by six independent param- 
eters. In the current fit we vary these six parameters of 
the direct 27r-production amplitudes. 

In the fitting procedure described above, nine one-fold 
ir + 7r~p differential cross sections are computed with non- 
resonant variable parameters of the JM model obtained 
by employing unrestricted normal distributions around 
their start values with cr-parameters of 10%-20% of their 
start values (see Table [TTJ) . In this way we explore mostly 
the range of ~ 3cr around the start parameter values. 
These ranges for non-resonant parameter variations in 
the JM model are shown in Table |TT] as a percentage of 
their start values. 

In this fit we also vary the "f v pN* electrocouplings 
and the ttA and pp hadronic partial decay widths of the 
Pn(1440) and Z?i3(1520) resonances around their start 
values. The start values of the Pn(1440) and £>i 3 (1520) 
electrocouplings are determined by interpolating the re- 
sults from the analysis [|0] of the CLAS data on A^7r elec- 
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troproduction off protons into the range 0.4 GeV 2 < Q 2 
< 0.6 GeV 2 and extrapolating these results into the Q 2 - 
area from 0.25 GeV 2 to 0.4 GeV 2 . The electrocouplings 
of the Pn(1440) and Di 3 (1520) resonances are varied em- 
ploying normal distributions with a parameters equal to 
30% of their start values. There are no restrictions on 
the minimum or maximum trial electrocoupling values. 
The normal distributions allow us to explore mostly the 
area of ~ 3cr around their start values or 90% around 
electrocoupling start values. 

The tt + tt~p electroproduction channel has also some 
sensitivity to the 511(1535) state, which couples domi- 
nantly to the Ntt and Nrj final states. The 511(1535) 
electrocouplings were taken from the CLAS analysis of 
Ntt electroproduction [l5| and varied strictly inside the 
uncertainties reported in that paper. 

The 7rA and pp hadronic decay widths of the 
Pn(1440), Di 3 (1520), and 511(1535) resonances are var- 
ied around their start values taken from previous analy- 
ses of the CLAS double pion electroproduction data [46l — 
l48l | in ranges restricted by the total TV* decay widths 
and their uncertainties, as shown in Table UTT1 The total 
N* decay widths were obtained by summing the par- 
tial widths over all decay channels. Partial hadronic 
decay widths to all final states other than ttA and pp 
are computed as the products of RPP [53| values of 
N* total decay widths and branching fractions for de- 
cays to particular hadronic final states, that are taken 
from analyses [56| of hadroproduction experiments. We 
varied the ttA and pp hadronic decay widths of the 
Pn(1440), L>i 3 (1520), and 5n(1535) resonances simul- 
taneously with their masses, keeping the hadronic N* 
parameters independent of Q 2 . Accounting for the cor- 
relations of the N* electromagnetic and hadronic decay 
parameters in a combined variation is important for a 
credible extraction of the resonance parameters and, in 
particular, for the evaluation of their uncertainties. 

The here developed fit procedure allows us to de- 
termine the 7rA and pp hadronic decay widths of the 
Pn(1440), L»i 3 (1520), and S"n(1535). The resonance to- 
tal decay widths are affected by the variation of their tt A 
and pp partial decay widths. Instead, when varying reso- 
nance electrocuplings, the N* total decay widths of reso- 
nances are kept almost unchanged. Resonance hadronic 
decay amplitudes are independent of Q 2 , while resonance 
electrocouplings represent the functional dependencies on 
the photon virtualities. These distinctive features allow 
us to disentangle resonance electromagnetic and hadronic 
decay amplitudes in the fit of the CLAS data @ and 
to obtain the 7rA and pp hadronic decay widths of the 
Pu(1440), Pj 3 (1520), and 5 n (1535). 

To account for the contributions from the tails of 
higher mass resonances, electrocouplings of all N* states 
marked in Table U as "var" arc varied around their start 
values. For resonances with masses above 1.6 GeV, ex- 
cept the 5n(1650) state, the start values for electrocou- 
plings are taken from the results of previous analyses 
[46|,|47j of tt + tt~p electroproduction data [43| after the 



Q' 2 bins, 


W intervals, 


Number of fitted 


X 2 /d.p. intervals 


(GeV 2 ) 


(GeV) 


points 


for selected 








cross sections 


0.25-0.30 


1.44-1.56 


378 




0.30-0.35 


1.41-1.54 


378 


2.66-2.74 


0.35-0.40 


1.41-1.54 


378 




0.40-0.45 


1.39-1.51 


373 


1.87-2.04 


0.45-0.50 


1.39-1.49 


309 




0.50-0.55 


1.39-1.51 


360 


1.57-1.75 


0.55-0.60 


1.34-1.44 


279 





TABLE IV: Quality of the CLAS n + n'p fit within the frame- 
work of the JM model Q ■ Horizontal lines separate the three 
Q 2 -intervals, in which the fits were carried-out independently. 



extrapolation into the Q 2 range of the current analysis 
(0.25-0.6 GeV 2 ). The 5u(1650) state may have a more 
pronounced impact on the extraction of the Pn(1440) 
and Pi 3 (1520) resonance parameters, since the 5n(1535) 
and «Sii(1650) states are mixed in the dressed resonance 
propagator of the unitarized BW ansatz employed in the 
JM model. For this reason we prefer to have more accu- 
rate start values for the 5n(1650) electrocouplings. They 
are computed within the framework of the single quark 
transition model (SQTM) [58[ from electrocouplings of 
the Sii(1535) state that are available from analyses of 
Ntt and Nr] electroproduction data [l[ . Electrocouplings 
of the L>i 3 (1700) and P 33 (1700) resonances are kept fixed 
at their start values, since the masses of these resonances 
are far outside of the W area covered by the CLAS 
tt + tt~p electroproduction data @ W < 1.6 GeV. 

The 7rA and pp partial hadronic decay widths for 
N* states other than the Pn(1440), P i3 (1520), and 
5ii(1535) resonances are taken from previous analyses 
of the CLAS double pion electroproduction data [4ftl48| . 
They are in a reasonable agreement with the values re- 
ported in the RPP [53]. The respective branching frac- 
tions are listed in Table fl] 

We fit the CLAS data @ consisting of nine differential 
cross sections of the ep — > e'p'TT + TT~ electroproduction 
reaction in all bins of W and Q 2 in the kinematic region: 
1.35 GeV < IV < 1.57 GeV and 0.25 GeV 2 < Q 2 < 0.6 
GeV 2 within the framework of the fit procedure described 
above. Three intervals of Q 2 separated in Table IIVI by 
horizontal solid lines are fit independently. The \ 2 /d.p. 
intervals, that correspond to equally good data descrip- 
tion within the error bars are shown in Table IIVI Their 
values demonstrate the quality of the CLAS tt + tt~p data 
description achieved in the fits. Examining the descrip- 
tion of the nine one- fold differential cross sections, we 
found that the y 2 j d.p. values are determined mostly by 
the random deviations of some experimental data points 
from the bunches of computed fit cross sections. There 
are no discrepancies in describing the shapes of the dif- 
ferential cross sections, which would manifest themselves 
systematically in neighboring bins of W and Q 2 . Typical 
examples for IV=1.51 GeV and neighboring Q 2 intervals 
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centered at 0.38 GeV 2 and 0.42 GeV 2 are shown in Fig.g] 
and Fig. respectively. 

The sets of computed differential cross sections with 
X 2 /d.p. values within the intervals given in the Table HVl 
offer the best data descriptions achievable within the 
framework of the JM model. The minimal values of 
X 2 /d.p. in each interval represent the global minima of 
these fits. We found that the increases of the \ 2 /d-P- 
values within the intervals listed in the Table IIVI change 
the computed differential cross sections, but still keep 
them inside the data uncertainties, offering equally good 
data descriptions in all these fits with different sets of JM 
model parameters. The mean values of the resonance pa- 
rameters from these sets and their dispersions (statistical 
r.m.s.) are determining the resonance parameters and 
their uncertainties, as it was described at the beginning 
of this Section. 

Since only statistical data uncertainties are used in 
the computation of the \ 2 /d.p. values listed in the Ta- 
bic IIVI we concluded that a reasonable data description 
was achieved. The x 2 1 d.p. values of our fits are com- 
parable with those obtained in the fit of the CLAS Nir 
electroproduction data published in II 511 . as well as with 
those obtained in the MAID analysis j62| . 



IV. EVALUATION OF ^ v pN* 
ELECTROCOUPLINGS AND RESONANCE 
HADRONIC DECAY PARAMETERS 

The extraction of resonance electrocouplings using the 
JM model relics on fitting resonant and non-resonant con- 
tributions to measured differential cross sections. There- 
fore, we first have to check the quality of the separa- 
tion between resonant and non-resonant contributions 
achieved in the data fit. 



A. Separation of resonant and non- resonant 
contributions to the n + n~p cross sections 

With the parameters determined from the fits, we now 
use the JM model to evaluate the contributions from res- 
onant and non-resonant parts to the cross sections. This 
is done by computing the nine differential cross sections 
without the resonant parts and with only the resonant 
parts for all trial differential cross sections of the JM 
model selected in the fit. In this way we also determine 
the ranges of the resonant and non-resonant contribu- 
tions to the cross sections as imposed by the uncertain- 
ties of the experimental data. These ranges account for 
the uncertainties of the non-resonant parameters listed in 
Table |TT]and for all correlations between resonance and 
non-resonant amplitudes. Therefore, we obtain reliable 
estimates for the uncertainties of the resonant and non- 
resonant contributions to the differential cross sections. 
Examples of the separated resonant and non-resonant 
contributions for two particular W and Q 2 bins are shown 



in Figs. 0] and [5] The interference between the resonant 
and non-resonant amplitudes is clearly seen in all angular 
distributions. The differences between the fitted ir + n~p 
differential cross sections and their non-resonant parts 
are larger than the resonant contributions for a majority 
of the data points in the angular distributions of the final 
hadrons, offering clear evidence for the interference be- 
tween resonant and non-resonant amplitudes. This inter- 
ference amplifies the N* contributions to all CM-angular 
distributions dcr/d(— cos(9i)) (i = 7r + ,7r~,j)) of the fi- 
nal hadrons, which improves the sensitivity of the fit to 
the resonance parameters. Furthermore, the shapes of 
the resonant and non-resonant differential cross sections 
are rather different, especially for the final-state angu- 
lar distributions. Substantial differences in the shapes of 
the resonant/non- resonant contributions and their inter- 
ference allow us to isolate the resonant contribution in 
a combined fit of all nine one-fold differential cross sec- 
tions despite the relatively small resonant contributions 
to the fully integrated cross sections. In our previous 
analysis [fj wc found resonant contributions to the fully 
integrated cross sections from 10% to 30%. 

The examples of Figs. @] and [5] demonstrate that the 
uncertainties of the resonance parts are comparable with 
those of the experimental data. This is further evidence 
for the credible separation between the resonant and non- 
resonant contributions. Any ambiguities in the evalu- 
ation of these two contributions would result in larger 
uncertainties for the resonant and/or non-resonant parts 
derived from the fit than the uncertainties of the original 
data. But this is not the case in the entire kinemati- 
cal area covered by the CLAS tt + tt~p electroproduction 
measurements of Ref. [f|. Therefore, these data provide 
enough constraints in order to determine the resonant 
contributions to all differential cross sections. 



B. iV* parameters from the fit of the tt + tv p 
electroproduction cross sections 

The procedure described in Section [TTT1 allows us to ex- 
tract the "f v pN* electrocouplings and their uncertainties. 
A special approach was developed for the evaluation of 
the Pi i (1440) electrocouplings. According to the anal- 
ysis [15| of the CLAS Nn electroproduction data, the 
A 1 / 2 electrocoupling of this resonance changes sign be- 
tween Q 2 of 0.40 GeV 2 and 0.45 GeV 2 . Here the A 1/2 
electrocoupling is close to zero and is an order of mag- 
nitude smaller than Si/%. The A l / 2 variations computed 
as a percentage of the start value, which is close to zero, 
become too small. For realistic uncertainty estimates we 
varied A l / 2 for 0.4 GeV 2 < Q 2 < 0.5 GeV 2 in the ranges 
shown in Table [V] By varying A 1 / 2 inside these widened 
ranges, we scanned the trial values comparable with those 
°f Si /2 , as they were obtained in the analysis [l5[ of the 
CLAS Nit electroproduction data. We refit the CLAS 
data of Ref. [|| on n + n~p electroproduction by vary- 
ing Ai/2, as described above, keeping the variation of all 
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FIG. 4: (color online) Resonant (blue triangles) and non-resonant (green open circles) contributions to the differential cross 
sections (red lines) obtained from the CLAS data 0], fit within the framework of the JM model at W— 1.51 GeV, Q 2 =0.38 
GeV 2 . The solid blue and dotted-dashed green lines stand for the resonant and non-resonant contributions, respectively, which 
were computed for minimal x 2 /d.p. achieved in the data fit. The points for resonant and non-resonant contributions are shifted 
on each panel for better visibility. Dashed lines show selected fits. 
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FIG. 5: (color online) The same as in Fig. H at W= 1.51 GeV, Q 2 =0.43 GeV 2 . 



other resonant and non-resonant parameters as described 
m Section Hill 

In order to compare our results on the Pn(1440) and 
Di3(1520) electrocouplings in the ir + ir~p electroproduc- 
tion channel with their values from the analysis of Ntt 



electroproduction [lj|, we have to use in both of the 
exclusive electroproduction channels common branching 
fractions for decays of these resonances to the Ntt and 
Ntttt final states. According to the RPP [53|, the sum 
of the branching fractions into the Ntt and Ntttt final 
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Q 2 bins, 
(GeV 2 ) 


Ranges covered 
in variations of Ax/2 
electrocoupling of the 
Pn(1440) resonance, 

(10" 3 GeV- 1/2 ) 


0.40-0.45 
0.45-0.50 


-20 - +5 
-5 - 20 



TABLE V: Variation range for the Pn(1440) A 1/2 electrocou- 
pling at photon virtualities where this electrocoupling changes 
sign. 



Q 2 , 


Ax/2, 


Sl/2, 


(GeV 2 ) 


(HP 3 GeV" 1/2 ) 


(10 -3 GeV~ 1/2 ) 


0.28 


-36.3 ± 17.0 


37.8 ±6.3 


0.33 


-21.5 ±7.7 


35.8 ±4.7 


0.38 


-10.9 ±5.9 


32.7 ±5.3 


0.43 


-8.8 ±5.4 


33.8 ±3.8 


0.48 


3.7 ±5.9 


37.4 ±5.9 


0.53 


15.8 ±9.1 


26.8 ±3.5 


0.58 


17.2 ± 10.4 


22.2 ±6.8 



TABLE VI: Electrocouplings of the Pn(1440) resonance de- 
termined from this analysis of it + tv~p electroproduction off 
protons within the framework of the JM model 0. 




Q (GeV 1 ) 




0.4 0.5 
Q 1 (GeV 2 ) 



FIG. 6: (color online) Electrocouplings and full error bars 
of the Pn(1440) state. The results from the analysis of the 
CLAS tt + tv~p electroproduction data [(| are shown by trian- 
gles. Electrocou plin gs from the analysis of the Nn electro- 
production data [3| are shown by circles. 



states accounts for almost 100% of the total decay widths 
of the P n (1440) and Di 3 (1520) states. In our analy- 
sis the branching fractions for P 1:L (1440) and £>i3(1520) 
resonance decays to the Nir final state BF(Ntt) were 
taken from theprevious CLAS studies of Nir electropro- 
duction data [HI], since the Nir exclusive electroproduc- 
tion channels are most sensitive to contributions from 
the Pn(1440) and L>i3(1520) resonances. The unitarity 
condition allows us to estimate the branching fraction 
BF(Nirir) corr value as: 



BF(Nirir) 



1 - BF(Nir) 



(14) 



For these resonance decays to the Nirir final states 
it turns out that the estimated branching fractions 
BF(Nirn) corr from Eq. are slightly (< 10%) dif- 
ferent with respect to those obtained from the ir + ir^p fit 
(BF(Nirir)o). Therefore, we multiplied the irA and pp 
hadronic decay widths of Pn(1440) and P>i3(1520) from 
the ir + ir~p fit by the ratio B ^[^Z 7T ' >c \' rr in order to make 
them consistent with the unitarity condition of Eq. (|14[) . 

Consequently, the Pn(1440) and L>i 3 (1520) electro- 
couplings obtained in our analysis are multiplied by cor- 
rection factors 



Chd — 



' BF(Ntttt) 
BF{Nirir) cor , 



(15) 



in order to keep the resonant parts and the full computed 
differential tt + tt~p cross sections unchanged under the re- 
scaling of resonance hadronic decay parameters described 
above. 

The Pn(1440) and £»i 3 (1520) electrocouplings derived 
from the fit of ir + ir~p electroproduction data are pre- 



sented in Tables fVTl and IVII1 The masses and hadronic 
decay parameters are shown in Tables [Villi and IIXI The 
hadronic resonance parameters are taken as the average 
of their values obtained in independent fits in the three 
Q 2 intervals shown in Table IIV1 

In Figures [S] and [7] we compare the results of this anal- 
ysis with results of the analysis of CLAS Nir electro- 
production data [HI]. Furthermore, there are also the 
Pn(1440) and Z?i3(1520) electrocouplings available from 
the MAID partial wave analysis of Nir electroproduction 
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Q (GeV 1 ) 




Q (GeV 1 ) 

FIG. 7: (color online) Electrocouplings and full error bars of 
the Di3(1520) state. This analysis of the CLAS -K + ir~p elec- 
troproduction data 0] is shown by triangles. Electrocouplings 
from the analysis of the Ntv electroproduction data [l5[ are 
shown by circles. 




Q 2 (GeV 2 ) 



FIG. 8: (color online) Comparison between the final results 
on 5*1/2 electrocouplings for the Pn(1440) from this analysis 
of CLAS 7T + 7r _ p electroproduction data [f| (triangles) with 
the values obtained by employing a regular BW ansatz for 
the resonant amplitudes and with the contributions from the 
S'ii(1535) state turned off (squares). 





Al/2, 


51/2, 


-43/2, 


(GeV 2 ) 


(10~ 3 GeV _1/2 ) 


(10~ 3 GeV _1/2 


(10~ 3 ) GeV" 1/2 ) 


0.28 


-51.9 ±4.8 


-44.0 ±4.1 


71.5 ±6.8 


0.33 


-57.1 ±7.2 


-44.5 ±4.1 


74.2 ±6.8 


0.38 


-73.4 ± 7.0 


-36.5 ±4.1 


65.3 ± 7.3 


0.43 


-75.0 ±9.3 


-43.9 ±5.3 


74.1 ±6.1 


0.48 


-82.6 ± 15.9 


-35.7 ±4.2 


63.0 ±6.8 


0.53 


-65.2 ±8.9 


-36.6 ±4.5 


59.6 ±10.0 


0.58 


-62.2 ±9.3 


-31.0 ±4.6 


48.9 ± 7.4 



TABLE VII: Electrocouplings of the Di3(1520) resonance de- 
termined from this analysis of -K + ir~p electroproduction off 
protons within the framework of the JM model 0] ■ 



at Q 2 < 0.6 GeV 2 [60l - l62| . Our analysis confirms the sign 
change of the A 1 / 2 electrocoupling for the Pn(1440), first 
observed in the Nir channels. We found the zero-crossing 
to be between Q 2 =0A GeV 2 and 0.45 GeV 2 in agreement 
with the Nir analysis (l5| . The electrocouplings for the 
Pn(1440) and £>i3(1520) are also consistent within their 
uncertainties with the CLAS analysis of the Nir electro- 
production [i~5| with the exception of the Ai/ 2 amplitude 
of Di3(1520), where we see a difference in the range of 
0.4 GeV 2 < Q 2 < 0.5 GeV 2 , which is slightly larger than 
a standard deviation. In this range of Q 2 two different 
analyses of the Ntt exclusive channels the CLAS and 
the MAID [6(|, give also slightly different values of the 
A 1/2 electrocoupling for the L»i 3 (1520). The MAID re- 
sults are close to our extraction at 0.4 GeV 2 < Q 2 < 
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Parameter 


Analysis of the CLAS 
7r + 7r~p data 


RPP 


Breit-Wigner mass, MeV 
Breit-Wigner width, MeV 
7rA partial decay width, MeV 

ttA BF, 
p p partial decay width, MeV 
p p BF 


1458 ± 12 
363 ± 39 
142 ± 48 
23%-58% 
6.2 ± 4.1 
< 2.0% 


1420-1470 (« 1440) 
200-450 (« 300) 

20%-30% 

< 8.0% 



TABLE VIII: Hadronic parameters of the Pn(1440) resonance determined from the CLAS data [f| on tt + tv p electroproduction 
off protons within the framework of the JM model 0] and from the RPP [k| ■ 



Parameter 


Analysis of the CLAS 
7r + 7r~p data 


RPP 


Breit-Wigner mass, MeV 
Breit-Wigner width, MeV 
7rA partial decay width, MeV 
ty A BF 

p p partial decay width, MeV 
p p BF 


1521 ± 4 
127 ± 4 
35 ± 4 

24%-32% 
16 ± 5 

8.4%-17% 


1515-1525 (w 1520) 
100-125 (w 115) 

15%-25% 

15%-25% 



TABLE IX: Hadronic parameters of the Di3(1520) resonance determined from the CLAS data [|J on 7r + 7r p electroproduction 
off protons within the framework of the JM model [3] and from the RPP [k§] ■ 



0.5 GeV 2 , but start to deviate at larger photon virtuali- 
ties. Instead, for the other electrocouplings of Di 3 (1520), 
both the CLAS [H| and the MAID [60[ analyses of the 
Ntt channels are consistent with the results of this anal- 
ysis. 

The unitarization of the resonant amplitudes and the 
contribution from the <Sii(1535) state have a negligible 
impact on A 1 / 2 for Pn(1440) and on all Di 3 (1520) elec- 
trocouplings. However, both factors affect the S\/2 elec- 
trocouplings of Pn(1440). Figure [5] shows a compari- 
son of the P\i(1440) S1/2 electrocouplings obtained with 
the unitarizcd BW ansatz incorporating the Sii(1535) 
contributions versus the regular BW ansatz without the 
<Sii(1535). The implementation of the Su (1535) contri- 
butions imposes additional constraints on the ranges of 
the Pn(1440) S1/2 electrocouplings that correspond to 
the computed differential cross sections, which are con- 
sistent with the data, allowing us to improve the uncer- 
tainties in the extraction of this electrocoupling. 

Consistent results on the electrocouplings of the 
Pn(1440) and Pi3(1520) resonances, which are available 
for the first time from independent analyses of the ma- 
jor Ntt and tt + tt~p exclusive channels provide evidence 
for the reliable extraction of these fundamental quanti- 
ties from the experimental data. Furthermore, this agree- 
ment strongly supports the reaction models employed for 
the analyses of the Ntt 15 1 and tt + tt~p Q electroproduc- 
tion data. 

The Breit-Wigner masses and mean values for the to- 
tal decay widths of the Pn(1440) and Pi 3 (1520), derived 
from the it + n~p data and listed in Tables NTU\ and |IX] 
are in a good agreement with those reported in the RPP 
[53]] . Uncertainties in the total resonance decay widths 



were obtained by varying the it A and pp partial decay 
widths, while the decay width to the Nir final state re- 
mained fixed at the values taken from Ref. [l5| . In com- 
parison with the results in the RPP (53|, our analysis 
gives somewhat larger branching fractions to the irA fi- 
nal states for the hadronic decays of both the Pn(1440) 
and Pi3(1520) states, and a smaller hadronic decay to 
the pp final state. 

Both jupN* electrocouplings and resonance hadronic 
decay widths are obtained at the resonant point on the 
real energy axis W = M r , using the relativistic Brcit- 
Wigncr paramctrizations of the resonant amplitudes. All 
resonance parameters determined in our approach incor- 
porate all relevant contributions to the N* structure. 
The extraction from the data of the individual contri- 
butions of the quark core and the meson-baryon dressing 
is outside of our scope. Our results can be compared 
directly with those determined from other meson elec- 
troproduction channels, where the BW ansatz for reso- 
nant amplitudes is used. For comparisons with results 
obtained from models that use other ways to describe 
the resonant contributions, the resonant cross sections at 
the resonant points can be compared. 



V. THE IMPACT ON THE STUDIES OF 
RESONANCE STRUCTURE FROM THE CLAS 
RESULTS ON THE Pn(1440) AND L>i 3 (1520) 
ELECTROCOUPLINGS AND HADRONIC 
PARAMETERS. 

In this section we discuss the impact of the CLAS 
data on the Pn(1440) and Di3(1520) electrocouplings 
and hadronic parameters determined from the indepen- 
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dent analyses of Ntt and ir + ir~p electroproduction off 
protons on the contemporary understanding of the struc- 
tures these states at various distances. 



A. Pii(1440) resonance 

The first comprehensive Pn(1440) electrocoupling sets 
available from the CLAS data on Ntt and tt + tt~p electro- 
production provide access to the active degrees of free- 
dom in the Pn (1440) structure at various distances. Pre- 
vious studies [lj, [l5| already showed a strong sensitivity 
of the electrocouplings to assumptions about the active 
components contributing to the Pn(1440) structure. 

The Pn(1440) resonance is characterized by peculiar 
features: 

• a smaller mass (1440 MeV) than that expected 
by quark models with vector flavor-conserving 
quark-quark (qq) interactions, this is less than 
the mass of the first orbital excitation of the nu- 
clcon, the Sn (1535), leading to so-called wrong 
Pii(1440)/Sii(1535) mass ordering [J, and 

• a large hadronic decay width of w 300 MeV, al- 
though other resonances with masses below 1.6 
GeV have decay widths of less than 200 MeV, 

that make it difficult to describe the Pn(1440) as a pure 
three quark bound state. To the best of our knowledge, 
there are no models, incorporating only quark degrees of 
freedom, capable of describing simultaneously the mass 
ordering of the Pn(1440) and Sii(1535) resonances, the 
unusually large Pi i (1440) hadronic decay width, and 
the photo and electrocouplings. Therefore, different ap- 
proaches for the description of the Pn(1440) were de- 
veloped that take into account the contributions from 
degrees of freedom other than those of dressed quarks 
Hi HHIl|. The CLAS results on the Pn(1440) elec- 
trocouplings of [TH and this analysis allows for the first 
time to confront the model expectations for the Pn (1440) 
structure with the experimental data, in particular at Q 2 
< 0.6 GeV 2 , where consistent results from the indepen- 
dent analyses of Ntt and tt + tt~p electroproduction off 
protons have now become available. 

In Ref. [13 the Q 2 -dependence of the P n (1440) elec- 
trocouplings was obtained by assuming a hybrid nature of 
this resonance. In this approach, the Pn(1440) is treated 
as a resonance dominated by a single SU(6) configuration 
of three dressed quarks oscillating against an explicitly 
excited glue, denoted q 3 G. The Pn(1440) as a hybrid 
state has a spin-flavor part of the wave function orthog- 
onal to that of the nucleon, whereas the spin-flavor part 
of wave function for the Pn(1440), as a radially excited 
three quark state, is identical to that of the nucleon. This 
difference makes it possible to distinguish between the 
Pn(1440) as a regular q 3 or hybrid q 3 G state by study- 
ing the <5 2 -evolution of the Pn(1440) electrocouplings. 
A characteristic prediction of [63] for the hybrid origin of 



the Pn(1440) is the absence of the longitudinal coupling, 
i.e. 5*1/2 = 0. This is a consequence of negligible SU(6) 
configuration mixing in the non-relativistic approxima- 
tion adopted in [63| for the description of the Pn(1440) 
structure. The recent Lattice QCD (LQCD) studies of 
the N* spectrum [H, HH also confirmed the presence of 
the leading SU(6)-spin-flavor-configuration in the struc- 
ture of Pn(1440). The S 1/2 = values for Pn(1440) 
are clearly in disagreement with the results of our anal- 
ysis of the Tr + ir~p electroproduction data (see Fig. [5]), 
showing that at small photon virtualities S , 1 / 2 is larger 
than A1/2, or at least comparable at Q 2 close to zero. 
This allows us to rule out sizable hybrid contributions to 
the Pn(1440), confirming the conclusions from our pre- 
vious analysis of the Ntt final state [l5j]. The lack of a 
substantial contribution of a hybrid component to the 
lowest mass resonances in the Pi i partial wave was also 
confirmed in the LQCD studies j34{ . 

The previous analysis [l5[ of the Ntt CLAS data 
showed that the most satisfactory description of the 
Pn(1440) electrocouplings was achieved at Q 2 > 1.5 
GeV 2 within the framework of relativistic light-front 
quark models [l(l[lj|, as shown in Fig.[9j In these mod- 
els, the Pn(1440) electrocouplings are evaluated from the 
fully covariant, one-body electromagnetic transition cur- 
rent for point-like constituent quarks in the impulse ap- 
proximation, e mployin g light-front relativistic Hamilto- 
nian dynamics [64| - l66| . Resonance electrocouplings are 
computed in a frame where the "+" component of the 
virtual photon light-cone momenta is equal to zero. For 
this particular choice, the contribution from qq pairs into 
photon propagation are eliminated (l3| . In these models 
JTot HH the Pn(1440) is treated as the first radial exci- 
tation of three constituent quarks. The wave function of 
the ground state and Pn(1440) are described under the 
simplifying assumption, that they are unmixed single os- 
cillator basis states. 

A reasonable description of the CLAS data at Q 2 > 
1.5 GeV 2 is achieved in both models [IS Ell! as seen m 
Fig. [HI Accounting for the relativistic transition operator 
is important in order to describe the CLAS data over the 
entire range of photon virtualities, and especially at small 
Q 2 where the sign of Ai/ 2 changes [Tol Il3l|. although a 
quantitative description of the data at small Q 2 has not 
been achieved (see Fig. [9]). 

A similarly good description of the Pn(1440) electro- 
couplings at Q 2 > 1.5 GeV 2 was found using a valence 
quark model [6 71 ] based on the covariant spectator formal- 
ism [68|, |69[ . The results of this model are also shown in 
Fig. IH1 Incorporating three constituent quarks only, this 
model treats the N* electroexcitation as a virtual pho- 
ton interaction with a valence quark in the ground state, 
while the two other quarks are combined to a spectator 
diquark in both spin states and 1. The wave function 
of the ground state is parametrized and fit to the data 
on elastic nucleon form factors. The Pn(1440) structure 
is described as the first radial excitation of the nucleon 
ground state quarks. With no additional parameters, 
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the wave function of the Pn(1440) was calculated em- 
ploying the orthogonality condition between the ground 
state and the Pn(1440) wave functions. The electrocou- 
plings of Pn(1440) are calculated for the valence quark 
transition between the nucleon ground state and its first 
radial excitation. As shown in Fig. [9j this model [67[ de- 
scribes the trend of the CLAS data for Q 2 > 1.5 GeV 2 . 

Different approaches to describe the Pn(1440) res- 
onance using structureless constituent quarks of non- 
running masses - two relativistic light front quark mod- 
els [Io|, [H| and a valence quark model [63] ™ are capa- 
ble of reproducing the major features of the CLAS data 
on P n (1440) electrocouplings only for Q 2 > 1.5 GcV 2 . 
These approaches fail to describe the electrocouplings for 
Q 2 < 1.0 GeV 2 . On the other hand, consistent results on 
Pn(1440) electrocouplings from the independent analy- 
ses of Nir and tt + it~p clcctroproduction off protons at 
Q 2 < 0.6 GeV 2 emphasize the credibility of the experi- 
mental data and the inability of describing the Pn(1440) 
structure by only accounting for quark degrees of freedom 
without meson-baryon dressing. 

The global analysis of Nir, photo-, electro-, and 
hadroproduction within the framework of the EBAC- 
DCC coupled-channel approach 0, OJ, HI, [ZQrlzl re- 
vealed this additional component in the N* structure 
that is most relevant at Q 2 < 1.0 GeV 2 . The general 
unitarity condition for full meson electroproduction am- 
plitudes requires contributions from non-resonant meson 
electroproduction and hadronic scattering amplitudes to 
the 'jvpN* vertex, as depicted in Fig. [TU] This con- 
tribution incorporates all possible intermediate meson- 
baryon states and subsequent meson-baryon scattering 
processes that eventually result in the N* formation. As 
was mentioned in Section QV] the ^ v pN* electrocouplings 
extracted in our analysis of tt + tt~p electroproduction 
data, as well as in the previous analysis of the Ntt data 
[l5j account for all contributions to the resonance struc- 
ture, including for those from meson-baryon dressing. In- 
stead, most of quark models, including aforementioned 
fiol fl3l loTj . account only for the quark-core component in 
N* photo-/electroexcitation, while meson-baryon dress- 
ing, being a part of the reaction mechanism, is completely 
outside of their scope. Therefore, the current and the 
previous CLAS results [lj| on ^ v pN* electrocouplings 
should be different with respect to those from the quark 
model expectations in Q 2 -arcas where the contributions 
from meson-baryon dressing are substantial. 

In Ref. [ill the meson-baryon dressing amplitudes are 
calculated from diagrams describing non-resonant mech- 
anisms in Ntt photo-, electro-, and hadroproduction 
channels with hadronic parameters fit to the world me- 
son hadroproduction data. The estimates of the absolute 
values for meson-baryon dressing contributions to the 
Pn(1440) electrocouplings showed that meson-baryon 
dressing amplitudes are maximal at Q 2 < 1.0 GeV 2 . 

Therefore the meson-baryon dressing contributions 
could be responsible for the differences observed in this 



area between the CLAS data on the N* electrocouplings 
and the quark model expectations of Refs. [ill EH l67| . 
The meson-baryon dressing decreases with Q 2 and in the 
region Q 2 > 1.0 GeV 2 we have a gradual transition to the 
dominance of quark degrees of freedom, as indicated by 
the better description of the Pn(1440) electrocouplings 
within the framework of quark models. 

Evaluations for bare-quark-core contributions to the 
dressed Pn(1440) electrocouplings have recently be ob- 
tained within the framework of QCD-based Dyson- 
Schwingcr Equations (DSEQCD) gl|. The DSEQCD 
results are close to the bare-quark-core contributions to 
the Pn(1440) electrocouplings extracted from the ex- 
perimental data within the framework of the EBAC- 
DCC coupled-channel approach [23| and far from the 
dressed Pn(1440) electrocouplings obtained in our anal- 
ysis. This DSEQCD result offers further evidence for 
substantial contributions from meson-baryon dressing to 
the Pn(1440) electrocouplings. 

We conclude that the structure of Pn(1440) is deter- 
mined by the combined contributions from an internal 
quark core of constituent quarks in the first radial exci- 
tation and an external meson-baryon dressing, the latter 
being most relevant at small Q 2 . 

Lattice QCD is making progress in calculating ^ V NN* 
transition form factors from first QCD principles. Ex- 
ploratory calculations of the 7„pPn(1440) transition 
form factors P 1 P 2 11 (Q 2 ) were carried out within LQCD 
in the quenched approximation using a pion mass of 
500 MeV and a simplified set of projection operators 
[2!||3(|. Despite the aforementioned simplifications, the 
first LQCD evaluations are in approximate agreement 
with the CLAS F[} X (Q 2 ) form factor results at Q 2 > 
1.5 GeV 2 . Improved LQCD calculations with dynami- 
cal quarks and a smaller pion mass of 380 MeV [3l], HH 
provide for the first time a reasonable description of the 
CLAS data for Q 2 < 1.0 GeV 2 . These studies demon- 
strated promising opportunities to relate the results on 
the Q 2 -evolution of the ^y v pN* electrocouplings directly 
to QCD. 

Together with the recent LQCD [HH3 and DSEQCD 
[Hi results on the N* spectrum and the first DSEQCD 
results on bare 7„pPii(1440) electrocouplings [4l|, this is 
an important step toward to our understanding of baryon 
formation and the nature of confinement from the first 
principles of QCD. 

The 7rA and pp Pn(1440) hadronic decay widths ob- 
tained in our analysis are listed in Table IVIIII and com- 
pared with RPP [53| values. The quark models, which 
describe N* hadronic decays through flux tube breaking 
[zH IzH, provide a good description of our ttA and pp 
decay widths of 100 and 2.5 MeV, respectively. However, 
the flux-tube breaking mechanism is unable to describe 
the Pn(1400) Ntt hadronic decays. It predicts a width 
of approximately 400 MeV for this channel, which is in 
strong disagreement with measured values [53|. Other 
approaches, listed in Ref. [l6|, fail to describe the Ntttt 
hadronic decay widths. 
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FIG. 9: (color online) Photo- and electrocouplings of the Pn(1440) resonance in comparison with predictions of quark models. 
The photocouplings are taken from RPP (open square) and the CLAS data analysis [HsJ (open triangle). Other data points 
are the same as in Fig. [6] The results of relativistic light-front quark models 0, [l3|] are shown by solid and dashed lines, 
respectively. Results of the covariant valence quark spectator diquark model are shown by the dashed dotted line. 
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FIG. 10: Processes that contribute to meson-baryon dressing 
in the N* photo/electroexcitations within the EBAC-DCC 
approach. 



To date the quantitative description of Pn(1440) 
hadronic decays remains a challenging problem. Con- 
sistent accounting for the meson-baryon dressing contri- 
butions to both electromagnetic and hadronic Pn(1440) 
vertices is important to g ain insight into the structure of 
this excited state l^lHlIflrTrill. 



B. Di3(1520) resonance 

Analyses of the previous CLAS data [HI on Z?i3(1520) 
electrocouplings determined from Ntt exclusive electro- 
production channels, showed that their satisfactory de- 
scription for Q 2 > 1.5 GeV 2 can be achieved within the 



framework of the hypercentral constituent quark model 
(hCQM) [73 (solid lines in Fig. \T$. In this model the 
central confinement potential is parametrized by a sum of 
Coulomb and linear terms expressed in collective hyper- 
coordinates for the three constituent quark system. The 
use of hypercoordinates effectively accounts for three- 
body effects in quark interactions. The remaining in- 
teraction between pairs of quarks is parametrized by a 
superposition of spin- and isospin-dependent hypcrfinc 
terms. Wave functions for resonances are obtained by 
diagonalizing the Hamiltonian in a non-rclativistic ap- 
proximation. Three parameters of the hQCM were fit 
to data on the baryon spectrum. With these parame- 
ters electrocouplings of all well-established excited nu- 
cleon states were evaluated for Q 2 < 5.0 GeV 2 , employ- 
ing non-rclativistic electromagnetic transition operators. 
The results obtained for the Z?i3(1520) state are shown 
in Fig. HU 

The hCQM model can describe the data trends for 
the £>i3(1520) electrocouplings for Q 2 > 1.5 GeV 2 . Pro- 
nounced differences between hCQM expectations and the 
Z?i3(1520) electrocouplings at smaller Q 2 offer an in- 
dication for contributions from active degrees of free- 
dom other than a core of three dressed quarks to this 
state's electrocouplings. Meson-baryon dressing ampli- 
tudes are likely contributors in this area of photon vir- 
tualitics. According to the EBAC analysis [17j . for the 
-Di3(1520) state they are maximal for small photon vir- 
tualities and decrease with Q 2 . The EBAC analysis also 
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FIG. 11: (color online) Description of the CLAS data on 
Di3(1520) electrocouplings within the hQCM (solid lines) 
f?l\ . The symbols associated with the experimental data are 
the same as in Fig. [9] 



suggests negligible meson-baryon dressing or dominant 
quark core contributions to the A 1 / 2 electrocoupling of 
the Di3(1520) state for Q 2 > 1.5 GeV 2 , where we already 
have the CLAS results on this electrocoupling [l5| . These 
results thus offer access to quark degrees of freedom in 
the structure of -Di3(1520) and open up new prospects to 
study the dynamical dressed quark mass, structure, and 
their strong interactions, that are responsible for the N* 
formation. These studies arc of particular importance 
to understand the nature of confinement in the baryon 
sector based on QCD [H M, ED, HI MM, III ■ 

Our analysis confirms a rapid helicity switch from the 
dominance of the A3/2 electrocoupling at the photon 
point to a comparable contribution from both transverse 
electrocouplings at Q 2 w 0.5 GeV 2 as already observed 
in [ll| . This is shown in Fig. [T2] in terms of the helicity 
asymmetry, defined as 



Ahel — 



A 2 



A 2 

^3/2 



A 2 



A 2 ' 

^3/2 



(16) 



This particular feature is expected for the contributions 
from the core of three constituent quarks in the first or- 
bital nuclcon excitation L = 1 [79]. It suggests a sig- 
nificant contribution from the core of three constituent 
quarks to the transverse Z?i3(1520) electrocouplings even 
at small Q 2 . 

We conclude that the Q 2 -evolution of the L>i3(1520) 
electrocouplings is consistent with contributions of both 
an external meson-baryon cloud and an internal core of 
three constituent quarks. 

The A^7r7r hadronic decay widths of £>i3(1520), ob- 
tained from the CLAS tt + tt~p data, are given in the Ta- 
ble ILXl They are close to the RPP [53J values. However, 
the branching fraction to the 7rA final state is larger and 
to the pp final state is smaller than the RPP average. Our 
results on the A^7T7r hadronic couplings of the Di3(1520) 
are well described by the flux tube breaking model of 
[zH lzE|: as wen as within the framework of the 3 Pq pair 
creation model of [Sfj]. The success of these models in- 
dicates a substantial role of the quark core in hadronic 
decays of the Di3(1520) state. 



VI. SUMMARY AND OUTLOOK 



The analysis of a large body of differential and fully 
integrated cross sections for the process "f v p — > n + Tr~p 
carried out within the framework of the JM model in our 
previous paper 0] allowed us to establish all essential 
mechanisms that contribute to this exclusive channel for 
1.3 GeV < W < 1.6 GeV and 0.2 GeV 2 < Q 2 < 0.6 GeV 2 . 
In this paper we use the reaction model of Ref. Q in or- 
der to determine the resonant and non-resonant contri- 
butions to the n + Tr~p differential cross sections measured 
with the CLAS detector @ and eventually to obtain the 
^1/2, S1/2, anc l A3/2 electrocouplings, as well as the 7rA 
and pp decay widths for the Pn(1440) and L>i 3 (1520) 
excited proton states from this data fit. 
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FIG. 12: (color online) Helicity asymmetry Ahei of Eq. (|16[) 
for the transverse electrocouplings of the Di3(1520) from this 
analysis of n + ir~p (filled triangles) and Ntt (open circles) ex- 
clusive electroproduction channels measured with CLAS. The 
curve represents the non-relativistic quark model calculations 



The good description of all measured observables 
makes it possible to evaluate the resonant contributions 
to the cross sections, which are needed to extract the 
j v pN* electrocouplings. Resonance electrocouplings and 
hadronic decay parameters were determined using a uni- 
tarized BW ansatz for the resonant amplitudes, that 
takes into account interactions between the same and 
different excited states in the dressed resonance propa- 
gators. 

For the first time electrocouplings of the Pn(1440) 
and Z?i3(1520) resonances were obtained from the anal- 
ysis of the 7r + 7r _ p exclusive electroproduction at photon 
virtualitics 0.25 GcV 2 < Q 2 < 0.60 GcV 2 . They are 
in reasonable agreement with the electrocoupling val- 
ues obtained in the previous CLAS analysis of the ex- 
clusive Ntt electroproduction channels |15| . Single and 
charged double pion electroproduction channels repre- 
sent two major contributors to meson electroproduction 
in the kinematical region covered by our measurements. 
The non-resonant mechanisms in these exclusive chan- 
nels are completely different. A successful description of 
a large body of observables in single and charged double 
pion electroproduction channels with consistent values of 
the Pn(1440) and Pi3(1520) electrocouplings confirms 
that the "f v pN* electrocouplings can be reliably deter- 
mined in independent analyses of these electroproduction 
channels. A good agreement between electrocouplings of 
the aforementioned excited states obtained from Ntt and 
tt + tt~p electroproduction also demonstrate that the reac- 



tion models developed to describe these exclusive chan- 
nels @, HH can be used to determine electrocouplings of 
most of the excited proton states that decay preferen- 
tially into either Ntt or Nirir final states. 

The 7rA and pp partial hadronic decay widths of the 
Pn(1440) and Di3(1520) states were also determined. 
They are close to the results of experiments with hadronic 
probes reported in the RPP [53| . 

The comparison of the CLAS data on Pn(1440) and 
£>i3(1520) electrocouplings with the expectations of con- 
ceptually different quark models [l(| EH an d , com- 
plemented by evaluations of meson-baryon dressing [ItJ 
allowed us to shed light on the relevant components in 
the structure of these states and to explore their evolu- 
tion with photon virtualitics. The electrocouplings reveal 
two major contributions, from an internal core of three 
dressed quarks and an external meson-baryon cloud. 

The structure of the Pn(1440) is consistent with a 
combined contribution from constituent quarks in the 
first radial excitation and from a meson-baryon cloud 
with pronounced contributions at Q 2 < 0.6 GcV 2 . The 
electrocouplings of the Pi3(1520) contain a combined 
contribution of three constituent quarks in the first or- 
bital nucleon excitation L = 1 and an external meson- 
baryon cloud, which becomes negligible for the A 1 / 2 elec- 
trocoupling of L»i 3 (1520) for Q 2 > 1.5 GcV 2 . The CLAS 
results on the Ai/ 2 electrocoupling of the Pi3(1520) res- 
onance in this kinematical region provide direct access 
to quark degrees of freedom in the structure of this ex- 
cited state. A physical interpretation of these results is 
of particular interest for those who are stud ying hadron 
structure starting from the QCD Lagrangian [28[ [H, l3lT - 

11, 51 til. 

The resonance electrocouplings and hadronic decay pa- 
rameters presented in this paper were obtained using the 
unitarized BW parametrization of resonant amplitudes. 
Their values are determined at the resonant point on the 
real energy axis (W = M r ) and incorporate all com- 
bined relevant contributions (quark-core, meson-baryon 
dressing and so on). The N* parameters extracted from 
the data in this way can only be compared directly with 
those obtained from the same or other exclusive channel 
data fits that employ BW parametrizations of resonant 
amplitudes. It remains to be seen how these resonance 
parameters can be related to the values determined from 
singularities of the reaction amplitudes in the complex 
energy plane, as employed in coupled-channel analyses 
[23| ■ fn the Section III CI we described the JM model re- 
lations between the resonance parameters and the model 
independent observables: the resonant part of the cross 
sections and N* electromagnetic/hadronic decay widths. 
We can use the model independence of these observables 
and require agreement of the cross section and/or reso- 
nance decay widths at the resonant point evaluated in 
different models. These observables, in turn, could be 
related to the resonance parameters for any particular 
model. This should make it possible to establish relations 
between resonant parameters determined in our approach 
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and those from a global mult i- channel analyses. 
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Appendix A: Resonance hadronic decay amplitudes 
employed in the JM model 

The relationship within the JM model [3, H3] be- 
tween the N* hadronic decay amplitudes (Xf\Td ec \^R) 
and the energy-dependent N* partial hadronic decay 
widths T\, (W) is presented in this Appendix. Here Xr 
is the helicity of the 'N*', which decays into irA or pp 
final states with helicity A / . 

The N* hadronic decay amplitudes in Eq. ^ can be 
expanded in partial waves of total momentum J: 

(A / |T dec |Afi) = ^(A / |T/ ec |Afi}< I/ (cosr)e«* ) (17) 



where 8* and cf>* are the CM emission angle of the tt for 
the 7rA intermediate state and of the final 'p' for the p 
p intermediate states. The indexes p and v are defined 
in Eq. ([3]). Only a single term with J=J ri where J r is 
resonance spin, contributes to the expansion in Eq. (|17[) 
because of total angular momentum conservation. We 
can rotate the quantization axis adopted for the initial- 
state | A.r) wave function and re-evaluate the matrix el- 
ement (Aj|Td ec |Aii) in the frame with the quantization 
axis defined by the final 7r (the final 'p') momentum for 
7rA {pp) N* decays, respectively. The matrix element 
(A/|T(i ec |Aii) after rotation becomes: 



decay width: 



(A/|T dec |Ail) = J2( X f\ T dec\Jr v')d J ; vl {cOs6* 



(18) 



The superposition of the states \J r v') in Eq. (TT5)) . with 
spin J r and projection onto the final-state quantization 
axis v' , is the transformed wave function of the initial- 
state | Aii) after the aforementioned rotation of the initial- 
state quantization axis. Rotational invariance of the res- 
onance hadronic decay amplitudes requires that the op- 
erator Tdec should be an SU(2) (g> 0(3) - spin ® orbital 
momentum scalar. Only the term with u 1 = v in Eq. 
(|18|) (with v defined by Eq. (|3])) remains non-zero in the 
sum of Eq. (|18j) , as a consequence of the Wigner-Eckart 
theorem applied to matrix elements (Xf\Td ec \J r v') with 
the scalar Td ec operator. 

From comparisons between Eqs. (|17[) and (|18[) we can 
see that: 



(h\TL\^) = (X f \T dec \J r u), 



(19) 



The (Xf\T^ c \X R ) matrix elements in Eqs. (|3TT9l) are de- 
termined by the final-state helicity A/ only, and are in- 
dependent of N* helicitics A.r. 

Assuming real values for the matrix element 
(A f\T^ c \X r ) in Eqs. (j3TTO|) . we can relate it to the T Xf (W) 
partial hadronic decay width of the excited state N* to 
the final state of helicity A / . We employ general relations 
of quantum theory between the resonance decay ampli- 
tude (A/|Xdec|Aii) of Eq. (fTTjl. the two-body state phase 
space of resonance decay products d$2&, and the T\ f (W) 



1 1 

2M r 2J r H 



?E / l(A/|T dec |Aii)| 2 d<i> 2fc .(20) 

Ail. 



The factor in Eq. 



20|) reflects the spin-tensor 
normalization in the convention of the JM model Q. 
This normalization, and the expression for the 'S'-matrix 
adopted in the JM model 0, defines the final-state two- 
body phase space d<fr 2 b as: 



(/<!> 



2/j 



1 Pf 

4tt 2 4M r 

W 2 
Ef = — 



sin(6*)d6*d4>*, 
M)-M), 



(21) 



2W 



p f 



E2 f 



where Ef,pf are the energy and momentum modulus of 
one of the final hadrons / (/ is either pion or the final 
proton for the N* — > irA or p 'p' decays, respectively), 
Mf is its mass, while the index /' stands for the other 
hadron. All frame-dependent kinematic variables of the 
final hadrons are defined in the final hadron CM frame. 
Inserting Eqs. (|1 71191) into Eq. (|20l) and accounting for 
and the 'd-'function normalization, 



J d'* v (cos9*) ■ d'l u (cos 8*) sin 8* d9* = 



2J+ 1' 



(22) 



we obtain Eq. ((3]) for the (A/|T(j ec |Aii) amplitudes, apart 

from the factor y j^y- Note that at the resonant point, 
W=M r , this factor is equal to unity. However, in calcula- 
tions of resonant cross sections for Wy^M r , the two-body 
phase space becomes different than that at the resonant 

point. The factor y^pj m Eq. accounts for this dif- 
ference. It is needed to evaluate resonant cross sections 
and amplitudes for W values off the resonant point. 

The (j>l) and (pi) absolute three-momentum values of 
the final ir for the N* —> irA decay (i=l) or of the final 
p' for the N* — > pp decay (i=2) in Eq. are averaged 
over the running mass of the unstable hadron in the in- 
termediate state: 



(pi) = J dM} 1 - 



M i0 r i0 



(Mf-Mf o r + M^Tl 



(Pi 



dMf- 



A^ n r 



iO 1 iO 



■ ' n(Mf-Mf ) 2 +Mf T 2 



P r i{M(), (23) 



Pi(M?), 



where Mi is the current invariant mass of the final irp 
particles in the case of N* — > irA decay (i = 1) or the 
current invariant mass of final irir particles for the N* — > 
pp decay (i = 2); Mio are the central masses of either 
A (i = 1) or p (i = 2); and r i0 , are their total decay 
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widths. The running momenta of the stable particles 
from N* decays Pi(Mf) in Eq. (|23|) are evaluated as: 

y 2W K 1 

[{W 2 + m 2 - M 2 ) 2 - W 2 ™ 2 ,] 1 ' 2 

P2 = 2W ■ 

The values of these momenta at the resonance point 
Pi(M?) were obtained from Eq. (JUJ) at W = M r . 

The hadronic decay widths r^. (W) were taken from 
experiments with hadronic probes, as was described in 
Section lim Those decay widths were obtained in another 
representation of orbital angular momentum L and total 
final hadron spin S TLS f (W). The partial T\ f (W) decay 
widths can be transformed into this representation [8l| 
by: 



2/-I-1 

^^ r <L05A 1 -A 2 |J,,A 1 -A 2 ) 
(si\ 1 s 2 - A 2 |5Ai - \ 2 ) 



where s\, Ai, s 2 , A 2 are spins and helicities for first stable 
and second unstable particles in the intermediate states. 



